Dynamics of a deformable body in a fast flowing soap film 



O 

o 

(N 

in 



q 

o 



> ■ 
O ^ 

o^ ; 
o , 
\o . 

o ■ 

o . 

o ■ 

• I— I 

Oh. 

> : 



Sunghwan Jung^, Kathleen Mareck-'^, Michael Shelley^, and Jun Zhang^'-'^ 
^Applied Mathematics Laboratory, Courant Institute of Mathematical Sciences, 
New York Universtty, 251 Mercer Street, New York, New York 10012, USA 
'^Department of Physics, New York University, 4 Washington Place, New York, New York 10003, USA 

(Dated: February 2, 2008) 

We study the behavior of an elastic loop embedded in a flowing soap film. This deformable loop 
is wetted into the film and is held fixed at a single point against the oncoming flow. We interpret 
this system as a two-dimensional flexible body interacting in a two-dimensional flow. This coupled 
fluid-structure system shows bistability, with both stationary and oscillatory states. In its stationary 
state, the loop remains essentially motionless and its wake is a von Karman vortex street. In its 
oscillatory state, the loop sheds two vortex dipoles, or more complicated vortical structures, within 
each oscillation period. We flnd that the oscillation frequency of the loop is linearly proportional 
to the flow velocity, and that the measured Strouhal numbers can be separated based on wake 
structure. 

PACS numbers: 47.32.ck, 47.54.De, 47.20.-k 



The wake flow behind a rigid obstacle is a central ob- 
ject of study in fluid mechanics. When the oncoming flow 
velocity exceeds a threshold, vortices are shed behind the 
obstacle Q . A typical wake is composed of successive ed- 
dies of alternating sign - the "von Karman vortex street" 
- and is observed over a wide range of flow velocities and 
body shapes 0, . The frequency of vortex shedding (/) 
is determined by the flow velocity {V) and the object size 
(d), whose relation is captured by the near constancy of 
the Strouhal number, St = d f /V 0. 

The dynamics of a rigid object which moves freely in 
the direction perpendicular to the flow is of interest in 
many industrial and biological applications Lat- 
eral motion of an object can be induced by interaction 
with the flow and is often called the vortex-induced vi- 
bration (VIV). At low flow velocities, the body starts 
to oscillate sideways with small amplitude (less than 0.4 
times body diameter). Its associated wake structure is 
again a von Karman vortex street. However, further in- 
crease of flow velocity causes the obstacle to oscillate in 
phase with the vortex shedding, and as a result, a series 
of dipoles are shed instead 0, ill • 

Settling bodies or rising bubbles, where the balance of 
gravitational and drag forces set the velocity, also exhibit 
transitions as they interact with their wakes. For ex- 
ample, a slowly settling sedimenting sphere falls straight 
downwards but above a certain sedimentation velocity, 
the sphere's motion becomes periodic and its trajectory 
a spiral or zigzag ^ deformable object, such as a 

droplet or bubble, can behave similarly even as its shape 
now changes [llLll^ . 

Finally, studies have shown the instability (and bista- 
bility) of slender deformable bodies to lateral oscilla- 
tions in quasi-2D soap-fllm flows [l^, and of heavy de- 
formable sheets to lateral oscillations in fast 3D flows 
0,0,0,0]. In these cases, the system corresponds 
to the flapping of a flag in a stiff breeze. 

Flowing soap film provides a practical template upon 



which to study the dynamics of a nearly 2D flow 0, . 
The e xper imental setup has been introduced earlier 
\lA list lia |2flj . In this work, we introduce a deformable 
closed body into a fast flowing soap-fllm. Two thin ny- 
lon wires (0.3 mm in diameter) separate at a nozzle (0.5 
mm inner diameter) attached to the bottom of a reser- 
voir. The reservoir contains soapy water maintained at 
a flxed pressure head, thus fixing the flux. A stopcock 
regulates the flow rate through the nozzle. The nylon 
wires extend downwards to a collection box 2.4 m be- 
low. Driven by gravity, the soap film flows downwards. 
Owing to air drag, a terminal velocity is reached approx- 
imately 60 cm below the nozzle with a velocity profile 
near the center close to being uniform (velocity differ- 
ences are within 20 % of the mean, over 60 % of the span 
about the midline). From optical interference patterns, 
the film thickness is found to vary smoothly across the 
film by about 15 % of its average thickness. 

We use a thin rubber loop (0.2 mm thick) as the de- 
formable structure. The loop wets into the soap-film and 
is supported from its inner side against the flow. The 
loop is much thicker when compared to the film thickness 
(0.003 mm) that the fluid presumably does not penetrate 
over the loop. Six loops of different circumferences (5- 
7.5 cm) are used. We find that for regimes studied here, 
the loop appears to undergo only bending deformations, 
and not stretching or compression, as its length shows 
no measurable increase or decrease. Currently, we do 
not understand what balance of effects sets the enclosed 
area of the loop, which is an important constraint on the 
possible dynamics. However, we do find that, once ex- 
perimental conditions are fixed, and the loop is in a fixed 
state of dynamics, the enclosed area changes very little 
in time (e.g. ~ 5% for a 5 cm loop). However, between 
different states or conditions, the enclosed area can vary 
by factors of two or three. 

A laser Doppler velocimeter (LDV; Model LDP-100, 
TSI Inc.) is used to record the upstream velocity V . 
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FIG. 1: Flow structures behind a 5 cm loop at 2.2 m/s flow ve- 
locity. The coupled fluid-structure system shows bistability: 
(a) the stationary state; the loop remains essentially motion- 
less and its wake is a von Karman vortex street. The loop 
is deformed by the flow into a teardrop shape, (b) the os- 
cillatory state; the loop sheds two vortex dipoles within each 
oscillation period. 

Micron-sized particles (Ti02) are seeded into the flow for 
LDV measurements. Flow structures are visualized using 
interference patterns from monochromatic illumination 
(low-pressure sodium lamps operating at wavelength 585 
nm). The movies of the wake flow together with the loop 
are recorded using a high speed camera at 1000 frames 
per second. 

The interaction between the loop and the flow is quite 
complicated. In our experiments, we observe bistable 
states, one stationary and another oscillatory (see Fig.^ 
and b), that co-exist over a range of flow velocities. At 
least in the conditions considered here, we do not observe 
spontaneous transitions between these two states. How- 
ever, a transition from the stationary to the oscillatory 
state can be induced by externally perturbing the loop, 
or by abruptly changing the flow velocity. 

In the stationary state, the loop behaves as a rigid 
hoop and has a teardrop shape with higher curvature 
on the top than on the sides (Fig.^). A characteristic 




FIG. 2: (a) At 2.2 m/s flow velocity, eight snapshots of a 5 cm 
loop and its centroid are shown. The oscillation period of the 
oscillating loop is 52 ms. The snapshots of the flexible loop are 
presented as solid loops at 6 ms intervals and the trajectory 
of the centroid of the loop with 2 ms intervals. The loops are 
sequentially numbered, (b) Wake structure when the loop 
centroid is at the far left. The loop starts to shed a counter- 
clockwise vortex. The wake structure is schematically shown 
in (c). 

length scale (D) of the loop, its width in the film, is 1 cm. 
The flow velocity (V) varies from 1.5 to 2.5 m/s and the 
kinematic viscosity (v) of soap film is approximately 0.04 
cm^/s. The frequency of vortex shedding (/^) varies from 
20 to 50 Hz. Based on these characteristic numbers, we 
estimate the Reynolds number (Re) and Strouhal number 
{Sts) for the system to be 

VD Df 

i?e= 5,000, St, = ^'^0.2, (1) 

v V 

where the subscript s stands for vortex shedding, since 
the Strouhal number is calculated based on the vortical 
structure of the wake. Figure^ shows the deforming 
body and its vortical wake using a 5 cm loop and flow 
velocity of 2.2 m/s. As can be seen, vortices of alternating 
sign are successively produced. 

In the co-existing oscillatory state, shown in Fig. 
at the same parameters as above, the loop now oscillates 
periodically in the horizontal direction, and the vortical 
wake behind it is quite different. For low flow velocity, 
two dipole pairs are shed during each oscillation period. 
Such a wake structure is also observed behind oscillating 
cylinders and is referred to as the 2P mode 8] . 

Figure|2jL shows both the position of the loop at several 
time-points during one period of oscillation, and the path 
taken by its centroid. During the oscillation, the loop 
continuously changes its shape, and its centroid moves 
along a figure-eight trajectory (Fig.Et)- This figure-eight 
shape is due to the fact that the frequency of oscillation in 
the stream-wise direction is twice that in the transverse 
direction. This has been observed in the motions of a 
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FIG. 3: Frequency of an oscillating loop (fioop) versus rescaled 
velocity {V^/S /a\^). We test over 6 different loop lengths (5 
cm (A), 5.5 cm (v)i 6 cm (t>), 6.5 cm (O), 7 cm (□), and 
7.5 cm (O))- The open symbols are in the oscillatory state 
and the closed ones are in the stationary state. The frequency 
of the oscillating loop is linearly proportional to the rescaled 
velocity. 



flapping flag 10 and VIV systems Also, the loop 

oscillates in phase with that vortex shedding (Fig.l^b and 
c). When the loop is at far right (or left), a clockwise 
turning (or counter-clockwise) vortex is shed. 

By using loops of several different lengths, we find a lin- 
ear relation between the oscillation frequency of the loop 
if loop) and a rescaled velocity VVS/a^/L (see Fig.OJ, 
where S is the film thickness, a the thickness of the loop, 
and L the loop length. Our results from loops of differ- 
ent lengths and differing flow velocities all collapse onto 
a single line with slope of about 0.27. This offset of this 
afhne relation suggests a bifurcation to oscillation at a 
finite flow velocity; linear extrapolation to fioop — sug- 
gests a critical rescaled velocity of about 20, which is 
unfortunately beyond the reach of this experiment. 

To better understand the relation between oscillation 
frequency and flow velocity, we propose a simple model 
for the oscillations of an elongated loop with longitudi- 
nal length Li driven by a "lift force" in the direction 
perpendicular to the stream. The lift force (F) is taken 
as proportional to pV'^LiS where p is the density of fluid, 
V the fluid velocity. Hence, F = {1/2)ClpV'^LiS, where 
Cl is a lift coefficient. Typically, L/ is proportional to 
the loop circumference, L. At an angle 9 inclined to the 
flow stream, Cl is proportional to sin 6* l^^- For small 
0, sin 9 ~ x^rrz/Vam whcre (a;^„,24„ is the center of mass 
(centroid) location. Therefore, we approximate the lift 
force as F = mx^ « {l/2)pV'^LiSx^^/y^^, where x^^ 
is the acceleration in the transverse direction and m is 
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FIG. 4: Strouhal number of the loop (StL) versus the flow 
velocity (V) for different loop lengths (5 cm (A), 5.5 cm (v)i 
6 cm (>), 6.5 cm (O), 7 cm (□), and 7.5 cm (O))- Open 
symbols are 2P mode and closed ones are flag-like modes. The 
Strouhal numbers of the two modes are well separated, and 
the Strouhal numbers of 2P modes are close to 0.2 whereas 
those of flag-like modes are above 0.25. 



the total body mass. In this experiment, the mass of the 
(wetted) loop is much greater than that of the enclosed 
fluid. Hence, we assume that the total body mass (m) 
is proportional to p^La? where pL is the density of the 
loop. Also, the y-component of the centroid, y^^ , and the 
length Li are assumed to be proportional to the length of 
the loop if the body is elongated due to the flow. With 
the trivial solution for the a:;-component of centroid as 
x^^^ — Ce"^*, we obtain an expression for the oscillation 
frequency: lo = 2Trf^^^^ oc VVS/a^/L. This is consis- 
tent with our observations and underlies our rescaling of 
the data in Fig. El Put differently, this is simply the 
oscillation frequency of a hanging pendulum where the 
gravitational force is replaced by a drag force. 

As the flow velocity increases, a more complicated 
mode in the oscillatory state can be observed (left panel 
in Fig.0J. In this case, the loop sheds more than four vor- 
tices over a single period of oscillation. We refer to this 
wake structure as a flag-like mode. We use this terminol- 
ogy because the wake now resembles more that behind a 
flapping flag (see 0), and because the body itself looks 
elongated and "flag-like". This is because the enclosed 
area is now smaller in relation to than for the body 
in the 2P mode. To characterize the loop oscillation and 
the wake structure, we again define a Strouhal number, 
now using the oscillation frequency of the loop itself, or 
StL = Af^^ /V where A is the outer amplitude of oscil- 



4 




FIG. 5: Enclosed area (normalized by /A-k) and amplitude 
(normalized by L) of the oscillating loop for 6 cm loop. Area 
and amplitude abruptly change as the wake structure tran- 
sits from 2P modes (open symbols) to flag-like modes (closed 
symbols). As the velocity increases further, the amplitude 
decreases due to the large drag force. 



lation as indicated in Fig.0] The corresponding Strouhal 
numbers for 2P modes and for flag-like modes are shown 
in Fig.0] The 2P mode yields St^ approximately 0.2 and 
the flag-like modes yield values above 0.25. 

This separation in the respective Strouhal numbers is 
caused by a discontinuous change of the amplitude A 
from the 2P mode to the flag- like mode. FigureEl shows 
the simultaneous transitions of enclosed area [S) and am- 
plitude {A) of the oscillating loop. At high flow velocity, 
the streamwise extension of the loop increases due to the 
higher drag, and this presumably causes the observed 
relative decrease in enclosed area in the flag-like mode. 
Due to its now higher aspect ratio, the loop's amplitude 
increases (compare left and right panels in Fig.^J). The 
causes of this abrupt change in enclosed loop area and 
oscillation amplitude remains an open question. Follow- 
ing this abrupt change, the frequency of vortex shedding 
also increases presumably due to the higher aspect ra- 
tio. Unlike the 2P mode, the oscillation of the loop and 
vortex shedding are not in-phase, and the wake becomes 
more complicated. 

We have reported on the dynamics of a flexible body as 
it interacts with an impinging high-speed flow. We find 
that the loop can have co-existing stable states. In the 
oscillatory state, the loop oscillation frequency is linearly 
proportional to the flow velocity and inversely propor- 
tional to the square root of the loop length. Also, it has 
been observed in a system with a freely vibrating solid 



cylinder that oscillation frequency linearly depends on 
the upstream flow velocity 0, ll| ■ This relationship can 
be explained by a simple model based on the lift force 
of an inclined airfoil. With small loop lengths and low 
flow velocity, a 2P wake mode (with pairs of shed vor- 
tices 0,0) is observed and has constant Strouhal number 
~ 0.2. At longer loop lengths and higher flow velocity, a 
flag-like mode appears with higher Strouhal number (> 
0.25). This separation in Strouhal numbers is a result 
of abrupt changes in enclosed loop area and oscillation 
amplitude. 
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